Mixed partial-wave scattering with spin-orbit coupling 
and validity of pseudo-potentials 
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We present exact solutions of two-body problem for spin- 1/2 fermions with isotropic spin- 
orbit(SO) coupling and interacting with an arbitrary short-range potential. We find that in each 
partial-wave scattering channel, the parametrization of two-body wavefunction at short inter-particle 
distance depends on the scattering amplitudes of all channels. This reveals the mixed partial-wave 
scattering induced by SO couplings. By comparing with results from a square-well potential, we 
investigate the validity of original pseudo-potential models in the presence of SO coupling. We find 
the s-wave pseudo-potential provides a good approximation for low-energy solutions near s-wave 
resonances, given the length scale of SO coupling much longer than the potential range. However, 
near p-wave resonance the p-wave pseudo-potential gives low-energy solutions that are qualitatively 
different from exact ones, based on which we conclude that the p-wave model can not be applied to 
the fermion system if the SO coupling strength is larger or comparable to the Fermi momentum. 



I. INTRODUCTION 

Two-body problem takes the most fundamental place 
in the process of exploring and understanding many- 
body properties. In particular, two-body solutions de- 
termine the essential interaction parameter in the mi- 
croscopic many-body Hamiltonian. In the field of di- 
lute ultracold atoms, the two-body interaction is gener- 
ally formulated by the zero-range pseudo-potential, pro- 
vided that it produces the same asymptotic two-body 
wavefunction at length scale much shorter than the mean 
inter-particle distance but longer than the range of real- 
istic potential. The generalized pseudo-potentials for all 
partial- waves were first derived by Huang and Yang[l|, 
and then improved later by Stock et aZ[2|. So far the 
most popular pseudo-potential is in s-wave channel de- 
scribed by a single s-wave scattering length, which can be 
improved by including the energy-dependence in a self- 
consistent way [HQ. Another popular one is the p-wave 
pseudo-potential, described by the p-wave scattering vol- 
ume which generally has strong energy-dependencejj, 

In view of the great success when applying pseudo- 
potential models to the homogenous or trapped atomic 
gases, it is generally believed that this model will equally 
apply to other configurations, such as in the presence of 
spin-orbit (SO) coupling. Recently, by sophisticated ma- 
nipulations of laser field and magnetic field, the NIST 
group has successfully realized an optically synthesized 
magnetic field for ultracold neutral atoms @. As a re- 
sult, an effective SO coupling is generated in the system 
along one direction. Subsequently there are several the- 
oretical proposals to realize the symmetric Rashba SO 
coupling [7], and it is conceivable that an arbitrary type 
SO coupling could be achieved in future experiments. As 
usual, all existing theoretical studies about the SO cou- 
pled system are carried out in the framework of s-wave 
pseudo-potential, i.e., using the s-wave scattering length 
as that without SO coupling(see recent review Q). Based 
on this model, the most remarkable effect of symmetric 



SO coupling is to support a two-body bound state with 
an arbitrarily weak interaction, due to the modified low- 
energy density of stateQ. 

Although pseudo-potentials have been justified under 
confinement potentials fH, [3, it is not obvious that it 
is still robust under the single-particle potential as spe- 
cial as SO couplings. In the two-body scattering process, 
trapping potentials and SO couplings have the same ef- 
fect in mixing different partial- waves, either due to the 
trap anisotrot)v[lol fTl| . or due to the intermediate cou- 
pling with spin sector. However, unlike the trapping 
potentials, which generally contribute a trivial constant 
potential as two particles get close, the SO coupling in- 
trinsically affects the kinetic term and thus still mix all 
partial-waves for the two-body wavefunctions at short 
inter-particle distance. This non-trivial effect is expected 
to have important influence on the validity of original 
pseudo-potentials in the presence of SO coupling. For 
instance, an obvious deficiency of original s-wave pseudo- 
potential is that this model can predict arbitrarily deep 
bound state with the binding energy scaled in terms of 
the SO coupling strength 9]; however, under a square- 
well (attractive) interaction potential the true binding 
energy must be lower-bounded by the potential depth. 
Moreover, this discrepancy can not be amended by taking 
into account the energy-dependence of s-wave scattering 
length, as we shall show later in Section IV. 

In this paper, we make efforts to exactly solve the 
two-body problem with SO coupling for a general short- 
range interaction potential, without resorting to pseudo- 
potential models. For simplicity but without the loss of 
essence, we have chosen the isotropic SO coupling and 
studied in the subspace where only s-wave and p-wave 
scatterings are relevant. We show that the short-range 
parametrization of the wavefunction in each partial- wave 
channel will additionally rely on the scattering amplitude 
of another partial- wave channel, which refiects the mixed 
scattering between different orbital channels induced by 
SO coupling. The exact form of wavefunction obtained 



2 



above allows us to solve the two-body problem under 
a square-well interacting potential. By comparing with 
results from s-wave or p-wave pseudo-potentials, we ad- 
dress the validity of the latter in the presence of isotropic 
SO coupling. We find the s-wave pseudo-potential pro- 
vides a good approximation for the low-energy scattering 
state and bound state solutions near s-wave resonance, 
with the correction depending on the strength of SO cou- 
pling, the finite range of the potential and contributions 
from p-wave channel. However, near p-wave resonance, 
using the p-wave pseudo-potential alone will lead to re- 
sults that are qualitatively different from exact solutions 
from the square-well potential. We conclude that the p- 
wave pseudo-potential can not be applied to fermion sys- 
tem if the SO coupling strength is larger or comparable to 
the Fermi momentum. We shall address the underlying 
reasons for these results. 

The rest of the paper is organized as follows. In sec- 
tion II, we present the exact solution for two spin-1/2 
fermions under a general short-range interaction poten- 
tial and with isotropic SO coupling. In section III we 
reduce the exact solutions to the framework of original 
s-wave and p-wave pseudo-potentials. In section IV we 
present the numerical results for two-body problem un- 
der the square- well potential, from which we address the 
validity of s-wave and p-wave pseudo-potential models. 
We summarize the paper in the last section. 



II. TWO-BODY PROBLEM WITH ISOTROPIC 
SPIN-ORBIT COUPLING 

In this section we shall solve the two-body scattering 
problem for a special case of isotropic SO coupling. As- 
suming a general form of short-range interaction poten- 
tial (See Eq|^ , we obtain the wavefunction of scattering 
statefEa l33)) and analyze its long-range and short-range 
asymptotic behaviors. Particularly we show that its 
short-range behavior is parameterized by scattering am- 
plitudes in all relevant partial- wave channels (Eqs l41l^ . 
Finally we present the bound state solutions which can 
be deduced from scattering state solutions via EaH51 

We start from the single-particle Hamiltonian of 
spin— l/2(t,^) fermions with isotropic SO coupling, (we 
set the reduced Planck constant Ti ~ 1) 



Hi = 



2m 



m 2m ' 



(1) 



where k = ikx,ky,kz) and cr — (cr^;, CTj^, CTz) respectively 
denote the momentum operator and Pauli spin opera- 
tor; A is the strength of SO coupling. The single-particle 
eigen-state has two orthogonal branches as 
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with (/)k = arg{kx + iky), u^'' = \J h ^ ^-^id the cor- 
responding eigen-energy e'^'^ = (|k| ± A)^/(2to) as shown 
in Figll] Due to the isotropy of SO coupling, the total 



angular momentum j = 1-t-s (s = ^cr) is conserved by Hi, 
giving the highest rotation symmetry among all types of 
SO couplings. 




FIG. 1: (Color online) Single-particle spectrum, e^'' — (|k| ± 
A)^/(2m), with isotropic SO coupling. For given energy eo — 
kg /{2m), two magnitudes of momentum are available as 
and k2, with fci = A — fco, k2 ~ X + kg. fco < A in (a) and 
fco > A in (b). 

The two-particle Hamiltonian can be written as H2 = 
Hji + i/k, with Hji and respectively describing the 
center-of-mass motion with total momentum K = ki -|-k2 
and relative motion with momentum k = (k2 — ki)/2. 
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With isotropic SO coupling, the total angular momen- 
tum for two particles J = L -I- S is also conserved, with 
L = li -I- I2, S = Si + S2 respectively the total orbital 
angular momentum and total spin of particle 1 and 2. 
Aloreover, L can be decomposed as L = Lr -t- Lr,, with 
Lr (Lr) the angular momentum for the relative motion 
r = r2 — ri (center-of-mass R = (ri -I- r2)/2). In view of 
the symmetry of H2, in this paper we consider the scat- 
tering problem in the subspace of K = 0(Lr, = 0) and 
J = Lj. + S = 0. (The method presented below can be 
generalized to the case of non-zero K or J). 

For total K = 0, the scattered wavefunction only de- 
pends on the relative coordinate r, and is given by the 
Lippmann-Schwinger eg nation [l^ as 



(r|*k) = (r|* 



(0)\ 



rfr'(r|G(i?)|r')(r'|C/|*k). (5) 



where G{E) = ^_J^^_^_^g is the two-particle green func- 
tion, U is the interaction operator; [vl/lf^) is the incident 
two-particle state with relative momentum k, which can 
be either of the following three states 

l^-)) - |k(-),-k(-))(-e^^>'), (6) 
|$(,++^) = |k(+),-k(+))(-e-^^>'), (7) 
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in coordinate space they are (we set the volume V = 1 
for normalization) 
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with fc = |k|, k± = yk"^ + k"^, k± — k^ ^ iky. 

Furthermore, the subspace of J = Lr + S = can be 
spanned by two orthogonal components as (labeled by 
\Lr, S;mL,ms)) 



\J = Q)o = |00;00), 



(12) 



|J = 0)i = — [|11;-1,1) + |11;1,-1)-|11;0,0)^13) 

Here | J = 0)o is the spin-singlet combined with s-wave 
orbital channel, and | J = 0)i is the spin-triplet combined 
with p-wave orbital channel. Now any state projected to 
J = subspace can be written as 



(r|^'k),/=o 
with bases 



{n,\j = o)o = roo(i^.) 



Mr){^r\J = 0)o + Mr^rlJ = 0)(il,4) 
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and wavefunctions 

Mr) = J dilr o(J = 0|a,)(r|%.); 
V'i(r) = / dilr i{J = 0\^lr){r\^k)■ 



m 

(17) 
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Here Ctr denotes the azimuthal angle of relative coordi- 
nate r, and Yim the spherical harmonics with azimuthal 
quantum numbers {l^m). After projected, the eigen- 
states of i?2, i-O-, Eas. (|9ll0llip . are given by 

(r|<'-))j=o = V4^bo(fcr)(f],|J = 0)o + 

iji(kr){nr\J = Q)il (19) 

(r|4++^)j=o = V4^[jo(fcr)(l],|J = 0)o- 

iji{kr){Vlr\J (20) 

(r|<-+))j=o = 0. (21) 

with ji{x){l = 0, 1) the spherical Bessel function of /— th 
order. Particularly, Eql^T] shows that ( — h) channel is 
not involved in the subspace of J = 0. 

Due to the single-particle spectrum modified by 
isotropic SO coupling (see Fig[T]), the incident state of 
two particles with energy E = k'^/m can be an arbi- 
trary combination of plane- waves with two different mag- 
nitudes of momenta, |k2| = A + /c and |ki| = |A — fc|. For 
k < X, 

(r|<°V=o - a(r|<i>(;-)),;=o + /3(r|<-))j=o;(22) 
and for k > X, 



(r|<°)),=o - a(r|<-')^.o + /3(r|$Lr'),/=o,(23) 
which both result in {k2 = A + fc, ki = X — k) 



(r|*[f^)j=o = VA^{[ajo{k2r)+Pjo{kir)]{nr\J^0)o + t[aji{k2r)+f3ji{kir)]{nr\J^0)i}. (24) 

In view of the property of H2, we also project the interaction U (with range ro) to J = subspace as 

(r|t/|*k),/=o = V4^[Fo(r)(ri,|J = 0)o + Fi(r)(17,|J = 0)i], (r < ro) (25) 

here Fq^ Fi denote the scattering amplitude in s-wave(ir = 0) and p-wave(Lr = 1) channel. The Green function in 
Eql5]is calculated by inserting a complete set of intermediate states fEalT51 and EO)) . 



(r\G\r')r l^ A^K''')^^^-'^) , {r\<^\:^>){'^>r>W) , 



Here the prefactor 1/2 is to eliminate the double counting of inserted states. 
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Combining Egs. ([5l [24l [25l [26l) . we obtain the closed form of scattered wavefunction (for r > tq) in each partial- wave 
channel (see Eq fTll as 



■00/^4-^- = ajo{k2r) + CfcJno(fc2?') - ^0(^2?')] + Pjoihr) + Cfc^ [no(fcir) +ijo{kir)], 
V'i/(i\/4^) = aji{k2r) + CkAni{k2r) ^ iji{k2r)] + I3ji{kir) + CkdMkir) + ijiihr)] (27) 



where {q — ki or k2) 



C, = ^^^^(/o(<z)-*/i(g)), (28) 



pro 

/o(g) = m drr'^Fo{r)jo{qr), 

JQ 

fi{q) = m drr'^Fi{r)ii{qr), 
Jo 



(29) 
(30) 



and ni{x) {I — 0, 1) the spherical Neumann fmiction of 
l—th order. 

We further simplify the complex wavefunction (j27p by 
employing the time-reversal symmetry, i.e., [H2,T] = 
where T is the time-reversal operator. Therefore we 
choose the wavefunction to be the eigen-state for both 
H2 and T. Noting that T{nr\J = 0)o = (rj^|J = 0)o, 
T{^r\J = 0)1 = —{ftr\J = 0)1, the only way to achieve 
= e*^* is to assume 



— —aainSe , 



(31) 
(32) 



with 5 = 6/2. Then up to a prefactor %/47r cos 5e"', Eal27l 
is reduced to 



"00 = Q;[jo(fef) - tan(5no(fc2r)] + 
P[io{kir) -I- tanfeo(fci7')], 
^i/i — a[ii{k2r) — ia,n5ni{k2r)] + 

/3[ji(^i?') +tanfei(fcir)]. (33) 

To this end we have obtained the exact form of scat- 
tered wavefunction for a given short-range potential U (r) 
defined in Ea l25l Eg 1331 reveals a unique scattering prop- 
erty in the presence of isotropic SO coupling, i.e., the 
wavefunction in each partial-wave channel is character- 
ized by two different momenta(see also FigH]) with oppo- 
site phase shifts. Note that without SO coupling, A = 0, 
k2 = — fci = k, Eql33] reduces to the standard form of 
s-wave and p-wave scattered wavefunctions in free space. 

The scattered wavefunction fEa lMl) has the following 
asymptotic behaviors at long-range and short-range of 
inter-particle distances. As fcr — > 00, the long-range be- 
havior is (up to a prefactor ^/4^^e^^) 



ipo 



siTi{k2r + S) sin(fcir — (5) 
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a- 



sin(/c2r — 7r/2 -t- S) 
k2r 



kir 



(34) 



sin(fcir - 7r/2 



At short-range tq < r ^ we have (up to a prefactor 
\/47rcos(5e'*) 



/3 + ( /3 , tan 5 
k2 ki r ' 



(36) 



ak2 + (3ki ,a /3 .tanS 
3 "+(^-^)^- 

For simplicity, we consider the limit of zero-range po- 
tential, i.e., assuming Fi{r) = j^Fi{r 0) {i = 0,1) 
in Eql25l Further according to Eqs. (|29|30p we introduce 



_ m — , 



0), /i = 



TO rFi{r — > 0) 
4^ 3 ' 



(38) 



whic h gives / o(g) = /g, fi{q) =_qfi [q = k2 or ki] 
Eqs. (|28l31l32p then relate /q and to a, [3, S as 



/o 
»7i 



sin5e*'^( 



g^aki 



/3k2~ 

kl • 
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(39) 
(40) 



Thus the short-range behavior fEqs l36l57| can be ex- 
pressed in terms of /q, f ^ as (up to a prefactor 47rcot 5) 



*/i(fc? + fci)-/o(fc? + fci) 



k2- 



^k2) 



— tan (5 

" Jo) ~; 



(41) 



3(fc2 - kl 



/o(fc?+M), + ,7^^(42) 



These results show that with SO coupling, the short- 
range parametrization of the wavefunction in each 
partial-wave channel will additionally depend on scat- 
tering amplitude of another partial-wave channel. This 
directly reflects the spin-mediated mixed scattering be- 
tween different orbital (partial-wave) channels, as is one 
of the most dramatic features of SO coupled system. 

At the end of this section, we study the bound state 
solution with energy E — —K^/m < 0. The bound state 
is given by the poles of scattering amplitudes (/o, fi — 
00), which corresponds to the following transformation 
from the scattering state f2!| 



k — ;> iK, 6 —ic 



(43) 



kir 



Using Eq|43l the bound state wavefunction can be de- 
duced from Eq |331 its long-range and short-range behav- 
iors can be deduced from Eqs. (l34l35p and Eqs. (|36l37p 
respectively. 
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III. PSEUDO-POTENTIAL MODEL IN 
INDIVIDUAL PARTIAL- WAVE CHANNEL 

The pseudo-potential model formulated in a given 
partial-wave channel is based on two assumptions. First, 
the interaction only acts on this particular channel. Sec- 
ond, the short-range behavior of wavefunction in this 
channel is still determined by the same scattering param- 
eter as that in the absence of SO coupling. The second 
assumption is based on a general belief as follows. If the 
range of interacting potential (rp) is much shorter than 
any length scale in the system, as inter-particle distance 
approaches r ^ , all other potentials are negligible in 
this limit and the asymptotic behavior of two-body wave- 
function is unchanged. The validity of pseudo-potentials 
has been verified in trapped systems in Ref. 0, H, Q . In 
the following we reduce the exact solutions obtained in 
Section II to the framework of s-wave and p-save pseudo- 
potential models. 



A. s-wave pseudo-potential 

The s-wave pseudo-potential corresponds to assuming 
Fi = 0, fi = 0; by mapping the short-range behavior 
of tpQ (Eq l4T1) to I/?- — l/cs with as the s-wave scattering 
length in free space, we obtain the phase shift as 

ta.n6 = -as ; . (44) 

k 

For scattering state, at low energies, tan (5 = — OgA^/fc 
giving the effective ID coupling giD = 2asA^/m, which 
is supported by the modified low-energy density of 
state(DOS) by isotropic SO couphngs (see also Ref.[§]); 
at high energies. Eg Ml reduces to tan (5 = —kus as in 3D 
free space. 

The equation for bound state solution is obtained from 
Eg Ml via transformations as EoH51 

- — K = - K^, (45) 
as 

which reproduces the result obtained by s-wave T-matrix 
approach^. Eq|45] results in a bound state solution 
for arbitrarily weak interaction, which is a direct con- 
sequence of the effective ID DOS at low energies. 



B. p-wave pseudo-potential 

The p-wave pseudo-potential corresponds to Fq = 
Oj/o = Oi ^-iid ^ is determined by mapping the short- 
range behavior of ibi (EgMj) to r/3 — Up/r^, with Vp the 
p-wave scattering volume in free space. We obtain 

^ ^ A4-H6A2fc2+fc4 

tand = —Vp . (46) 



Without SO coupling (A = 0), it reproduces the original 
free space result as tan (5 = —Vpk^. 

For scattering state at low energies, tan (5 — — I'px 
again giving S{k = 0) = 7r/2; at high energies, it recovers 
the free space result. 

For bound state, by transformation as Eg 1431 we obtain 
from EqHBj that 

- — K = A^ - 6X^K^ + n\ (47) 

Vp 

We see that for arbitrarily weak p-wave interaction Vp — !• 
0~, Eql47l gives a shallow bound state as k = —VpX^. 



IV. SCATTERING UNDER A SQUARE- WELL 
POTENTIAL AND VALIDITY OF 
PSEUDO-POTENTIALS 

In this section we present the scattering state and 
bound state solutions under a square- well interaction po- 
tential. By comparing these solutions with those from 
individual s-wave and p-wave pseudo-potential model, 
we shall address the validity of pseudo-potentials in the 
presence of isotropic SO coupling. In Appendix A we 
show more details about partial-wave scattering under 
the square-well potential without SO coupling, and in 
Appendix B we derive the equations for two-body solu- 
tions with isotropic SO coupling. 



A. Results 

We consider a square- well potential with depth Vo(< 0) 
at inter-particle distance r < tq and with depth zero 
otherwise. The interaction strength is uniquely charac- 
terized by a dimensionless parameter as gro, with q = 
V— toVq. Without SO coupling, Eq lA4l shows that by in- 
creasing gro, a sequence of s-wave resonances (with phase 
shift 5s = 7r/2) occur at qro/i: = n + 1/2 and p-wave 
resonances ((5p = t^/'^) at qro/n — n + 1 (n = 0, 1, 2...). A 
bound state emerges whenever across a scattering reso- 
nance. 

Next we solve the two-body problem in the presence 
of isotropic SO coupling. Based on exact solutions in 
section II, the wavefunctions inside the potential (r < tq) 
in orbital s-wave and p-wave channels are 

V'o = jo{q2r) +tjo{qir), 
V-iA = Ji{q2r)+tji{qir); (48) 

with q2 = X + ^/m{E-Vo), gi = A - ^/i^^UE^^Vo) . 
Outside the potential (r > rg), the wavefunctions are 
given by Eql33|for the scattering state(i^ = k'^/m > 0), 
or by the transformed form (through Eg (43)) for bound 
state {E = -K^/m < 0). 

Using the continuity properties of i/'Oi ipi and their 
first-order derivatives at the boundary r = rg, one can 
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solve all the unknown parameters {t, a, /3, 6} for the 
scattering state and {t, a, (3, k} for the bound state. In 
Appendix |B] we present the equations for these solutions. 
Next we show numerical results for the scattering state 
and bound state in turn. 



1. Scattering state 

For given energy E = /m > Q, we obtain two phase 
shift solutions, 5i with /q » /i and 82 with /i ^ fo, 
analogous to s-wave and p-wave phase shifts without SO 
coupling. For fixed SO coupling Arg = 0.2, we show in 
FigHIa) the solution of Si near the first s-wave resonance 
and in Figl^Jb) the solution of 62 near the first p-wave 
resonance. Independently, we obtain Si from Eg 1441 using 
s-wave scattering length(as) with effective-range correc- 
tions (see Eq lA3[ / = 0), and S2 from Eq|46] using p-wave 
scattering length(ap = Vp/rg) with effective-range cor- 
rections fEa lA31 / = 1). In Figl2l these results are shown 
(by orange dashed lines) to compare with exact solutions 
(black circles). 




0.0 0.1 0.2 0.3 0.4 0.5 
kr„ 



0.0 0.1 0.2 0.3 0.4 0.5 
kr„ 



FIG. 2: (Color online) Phase shifts of scattering states in 
a square-well potential with isotropic SO coupling strength 
Aro — 0.2. Exact solutions are shown (by black circles) 
in comparison with results from pseudo-potential models 
with effective-range corrections(see Ea lA3|) (by orange dashed 
lines). (a)(5i near s-wave resonance with (from bottom to top) 
qro/n = 0.25(as < 0),0.5(as = oo),0.7(as > 0). (b)52 near p- 
wave resonance with (from bottom to top) qro/n = 0.9(ap < 
0),0.99(ap < 0),l(ap = oo),1.02(ap > 0). 

For the solution 5i near s-wave resonance, Figl2l^a) 
shows that it can be approximately fit by s-wave model 
within Zero < 0.5. Particularly ai k = Q, 5i = 7r/2 is 
consistent with the s-wave prediction (Eg HI)) due to the 
ID feature of the low-energy DOS. However, there is still 
a small deviation between these two solutions at finite fc, 
due to the interplay between SO coupling, p-wave con- 
tribution and the finite potential range. To investigate 
these effects in detail, we further study the modified ef- 
fective scattering length aef / in s-wave channel, which is 
defined by ^o{f) — >■ I/?' — l/aeff{k) at tq < r ^ 1/k. 
Practically aeff{k) can be extracted from the asymptotic 



wavefunction (|4ip by diagonalizing Eg lBIl in Appendix B. 
FigEI^a) shows how Og/ / (k) evolves with k for each given 
SO coupling, which can also be expressed in the form of 
effective-range correction. 



Qeffik) aeff 



reffk^. 



(49) 



One can see that with increased SO couplings, the ef- 
fective range r^fj almost stay unchanged while l/ce// 
become smaller indicating weaker interactions. The de- 
viations of 1/aeff from l/a^ directly manifest the effect 
of SO coupling and mixed scattering of s-wave channel 
with p-wave channel. Moreover the mixing can also been 
seen from the additional dependence of ^po on the p-wave 
scattering amplitude in EqUT] In Figl2Ub) we show the 
zero-energy value 1 /fle/ / as a function of Aro for several 
different potential depths. At Aro ^ 1; can be 

well fit by 



- + C(Aro)^ 

Qeff as 



(50) 



where the dimensionless parameter C only depends on 
the properties of the potential, or the actual interaction 
strengths in s-wave and p-wave channels. In FigUl C is 
shown as a function of r^/ag (together with ap/us) near 
the first s-wave resonance. In the weak interaction limit, 
|Vo| — ?> and as,ap — > 0^, C change linearly with ro/os, 
indicating aeff — oc — (Aro)^ in this limit. For the typ- 
ical parameter regime in the present experiment @, A is 
determined by the wavevector of the laser which is much 
smaller than the cutoff momentum of realistic potential. 
In this case, the condition Aro ^ 1 gives negligible cor- 
rection to aeff near s-wave resonances. 
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FIG. 3: (Color online) (a)Inverse of effective scattering 
lengths {aeff{k)) as functions of k for different SO couplings 
Aro ~ 0, 0.1, 0.2, 0.3(from top to bottom). The potential 
depth is qro/n = 0.25. (b) Inverse of zero-energy scatter- 
ing lengths (aeff) as functions of A at different qro/n — 
0.25, 0.5, 0.7(from bottom to top). The lines are fits to Eq|50l 
with C = -2.10, -1.00, -0.58(from bottom to top). 

For the solution S2 near p-wave resonance, however, 
it behaves qualitatively different from that obtained 
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FIG. 4: (Color online) C in Eg 1501 as functions of ro/as near 
and across the first s-wave resonance (gro < it) in the square- 
well potential. The red dashed line is the linear fit in the 
weak interaction limit as C = —1.073 + 0.275ro/as. Inset 
shows the ratio of p-wave to s-wave scattering length, where 
the gray (light) lines denote the s-wave resonance. 



entirely in the framework of p-wave pseudo-potential 
model, as showed by Figl2jb). Obviously, the exact solu- 
tion shows the initial value S2{k = 0) — or tt, depending 
on whether or not there is a two-body bound state(see 
next section); while the p-wave model always predicts 
S2{k = 0) = tt/2 according to EqHHl We have checked 
that in the limit of Arg ^ 1, the exact solution of 62 at 
k ^ l/rp essentially follows the free space result (given 
by tan — —Vpk^) with ^ or tt; while the p-wave 
model gives a narrow momentum window as < fc < A 
when S2{k) evolves from tt/2 to the exact result. This 
dramatic difference indicates that even near the p-wave 
resonance, the p-wave pseudo-potential alone can not be 
applied to the fermion system if A is larger or comparable 
to the Fermi momentum. We shall analyze the reason for 
the breakdown of p-wave model to scattering state solu- 
tions in the discussion section. 



2. Bound state 

The bound state solution E — —K^/m < is given by 
the transformed matrix equation = 0(see Appendix 
B). By setting k = in the matrix equation we determine 
the critical potential depth |Vo|c, which is responsible for 
the emergence of a new bound state, by 



(51) 



with (72 = X + qc,qi ~ X — qc and qc — \/ m|Vo|c. The so- 
lution of (7c is shown in FiglS] As A approaches zero, one 
branch of solution(solid line) is given by ji{qcro) = 
or Us — 0; the other branch(dashed line) is given by 
joiqci"o) ~ or Cp = 00. For the first branch, when 
increasing A the lowest solution will stay at qcV^ = or 
fls = 0~, while the other solutions increase resulting in 
deeper potential depths. For the second branch, when in- 
creasing A all solutions of will decrease, implying that 



weaker interaction is required to support the new bound 
state near p-wave resonance. In all, we see that only 
the lowest solution of the first branch is consistent with 
the prediction from s-wave pseudo-potential model (see 
Eg 145]). but none of the other solutions. The discrepan- 
cies here are attributed to the mixed scattering between 
s-wave and p-wave channels induced by the isotropic SO 
coupling. 
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FIG. 5: Critical potential depth ((/cro/vr) for the emergence 
of each new bound state as a function of isotropic SO cou- 
pling strength. When Aro — ^ 0, the solid lines approach 
qcTo/iT = 0, 1.430, 2.459, ... with as = 0; the dashed lines ap- 
proach qcTo/iT = 1, 2, ... with Op = 00 (See text). 



As shown in Figl6] with fixed Arg ~ 0.2, a sequence 
of bound states will develop when the potential depths 
increase above critical |Vo|c. For comparison, we also 
present the results from s-wave and p-wave pseudo- 
potential models, using the scattering length with or 
without energy-dependence. (For the bound state, the 
energy-dependent scattering length is determined from 
Eas. (|A2IA4p but with k replaced by zk. [l]) 

Figini^a) shows that the s-wave model using s-wave 
scattering length without (os) or with(as{E)) energy- 
dependence both give good approximations to low-energy 
solutions near s-wave resonance, but deviate a lot from 
exact solutions for deep bound states. In general, we find 
that using as{E) provides more accurate results than us- 
ing as in a large energy-range; particularly, in the limit of 
zero SO coupling, using as{E) will give the exact bound 
state solutions Q]. For fixed potential depth, the de- 
viation of s-wave results from exact solutions increases 
with the SO coupling strength, as shown in Figl71[al,a2). 
Moreover, FiglHlJa) and FigHshow that the s-wave mod- 
els using as{E) always predict deeper bound states than 
real solutions, which is consistent with Ea lSOl and also 
FiglHlthat the presence of SO coupling reduces the effec- 
tive interaction parameter for low-energy states. 

In Figlll^b), we further plot the relative deviations, 
Ak/k, as functions of SO coupling strengths at differ- 
ent potential depths. It shows that Ak/k increases 
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FIG. 6: (Color online) Bound state solutions (solid lines) in 
the square-well potential as functions of the potential depths, 
in comparison with results from s-wave(a) and p-wave(b) 
pseudo-potential models. The SO coupling is Aro = 0.2. In 
(a), we have used s-wave scattering length without (green cir- 
cles) or with (black triangles) energy-dependence. In (b), we 
use the energy-dependent p-wave scattering length(blue tri- 
angles) . 



FIG. 7: (Color online) (al,a2) Bound state solutions in a 
square-well potential as functions of isotropic SO coupling 
strengths at potential depths qro/ir = 0.25[(al), as < 0] and 
0.7[(a2), as > 0], in comparison with results using s-wave 
scattering length without (red dashed line) or with (green dot 
line) energy-dependence, (b) Relative deviations between ex- 
act solutions and results from s-wave pseudo-potential model, 
at the same potential depths as in (al,a2). For s-wave pseudo- 
potential, we use s-wave scattering length without ('x') or 
with ('-|-') energy-dependence. 



more rapidly for deep bound states than that for shallow 
ones. As also mentioned in the introduction, the s-wave 
model(even using the energy- dependent as{E)) is quite 
questionable when applied to deep molecules. As shown 
in FiglSJ the energy of the bound state is always lower 
bounded by the potential depth Vq, i.e., k < q. How- 
ever, the s-wave model will produce unphysically deep 
molecules with k > q. In this case, the s-wave model 
alone will not work and one must take into account the 
effect mixed scattering with p-wave channel due to SO 
couplings. 



In Fig[6l^b) we show the comparison with results from 
p-wave pseudo-potential model. According to the p-wave 
model (see Section IIIB), the bound state exists for an ar- 
bitrarily weak interaction in the presence of isotropic SO 
coupling. This is qualitatively different from the exact so- 
lution under the square-well potential, where each emer- 
gence of a new bound state requires a potential depth 
beyond the critical value (as shown by red lines in FiglSJ. 
In the limit of Arp ^ 1, the critical depths continuously 
approach qcTo/n — 1,2,... as in free space. The break- 
down of p-wave model to bound state solutions will be 
discussed in the following section. 
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FIG. 8: (Color online) Bound state solutions as functions 
of isotropic SO coupling strengths for given potential depth 
gro/jr = 0.4, in comparison with results using s-wave scat- 
tering length without (pink circles) or with(black triangles) 
energy-dependence. The orange dashed line denotes the lower 
bound as K = g. 



B. Discussion 

Through the last subsection, we have shown that 
the SO coupling has different effects on the validity of 
pseudo-potentials in s-wave and p-wave channels. In the 
limit of Aro ^ 1, the s-wave pseudo-potential model pro- 
vides good approximations to the low-energy scattering 
state and bound state solutions near s-wave resonances. 
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For example, it predicts correctly the initial phase shift 
as Si{k = 0) = 7r/2 for scattering state, and a bound 
state solution for arbitrarily weak attraction — > 0^. 
However, near p-wave resonances the p-wave pseudo- 
potential will produce qualitatively different results com- 
pared with exact solutions. For example, in the limit of 
Aro <SC 1, the exact solutions approach free space results, 
i.e., (52(fc = 0) =0 or tt, and each bound state emerges 
when Vp goes across a resonance at certain critical po- 
tential depth; on the contrary, the p-wave model predicts 
S2{k = 0) = 7r/2 and a bound state for any weak p-wave 
interaction Vp ^ 0~ . 

Here we analyze the reason why the s-wave pseudo- 
potential is approximately valid for SO coupled sys- 
tem while p-wave is not. This can be explained from 
the correspondence between the assumptions of pseudo- 
potential models and the resulted short-range behavior of 
wavefunctions. For s-wave pseudo-potential(/o ^0, fi = 
0), the resulted wavefunction does not show l/r^ singu- 
larity in p-wave channel, which is consistent with the 
assumption of zero scattering amplitude in Eq l25l How- 
ever, the p-wave pseudo-potential(/o = 0, /i 7^ 0) will 
induce an additional singularity in s-wave channel, i.e., 
Tpoir) 2A[^^^f^-|-iS£-^] as r ^- 0. This in turn requires 
that the interaction operator U also generate scattering 
amplitude in the s-wave channel, contradictory with the 
initial assumption that /o = in Eq l25l This paradox 
also implies that any weak /o will have dramatic interfer- 
ence with the p-wave sector and lead to qualitatively dif- 
ferent results from those using pseudo-potential entirely 
in p-wave channel. 

The results presented in this section tell us that the 
original pseudo-potentials formulated in the absence of 
SO coupling do not necessarily apply to the case with 
SO coupling. The pseudopotential model will certainly 
breakdown if the results obtained are inconsistent with 
initial assumptions of this model. In this case, an appro- 
priate interaction model should be constructed in order 
to rightly incorporate mixed scatterings between differ- 
ent partial- wave channels, which is out of the scope of 
this paper. 



partial-waves that is induced by the SO coupling. 

(2) Under the conditions that the length scale of 
SO coupling much longer than the range of the 
potential(Aro ^ 1) and near s-wave resonances, the s- 
wave pseudo-potential gives a good approximation to the 
low-energy solutions, with the correction depending on 
the strength of SO coupling, the finite range of the poten- 
tial and contributions from other coupled partial-waves. 

(3) Near the p-wave resonance, the p-wave pseudo- 
potential model gives low-energy solutions that are quali- 
tatively different from the exact ones from the square- well 
potential. The p-wave model alone can not be applied 
to the fermion system when the SO coupling strength is 
larger or comparable to the Fermi momentum. Its break- 
down is attributed to the inconsistent treatment between 
the assumption of the p-wave pseudo-potential and the 
resulted short-range singularities of wavefunction in s- 
wave channel. 

Although above results are obtained for the special 
type of isotropic SO coupling, they reveal the generic 
scattering properties modified by the coupling between 
spin and orbit. We therefore expect these results have 
strong implications to other systems with a general type 
of SO coupling. 
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Appendix A: Free-space scattering lengths under a 
square-well potential 

Far away from the range(ro) of the potential, the scat- 
tered wavefunction in ^— th partial-wave channel reads 



= ii{kr) - i&n5ini{kr) 



(Al) 



here 5i is the phase shift which give the scattering length 
(with energy-dependence) as 



V. SUMMARY 

In this paper, we exactly solve the two-body problem 
of spin- 1/2 fermions with isotropic SO coupling under 
a general short-range interaction potential, and investi- 
gate the validity of s-wave and p-wave pseudo-potentials 
formulated in the absence of SO couplings. Our main 
results are summarized as follows: 

(l)In the presence of isotropic SO coupling, the two- 
body scattered wavefunction exhibits exotic dependences 
on the momentum and phase shift (Eg 1551) . In each 
partial-wave channel the wavefunction at short inter- 
particle distance is parameterized by scattering ampli- 
tudes of all coupled scattering channels (Eqs l41l4^ . This 
feature reveals the mixed scattering between different 



ai{k) 



1 tan 5i 



(A2) 



In the limit of kr^ <C 1, the effective-range expansion 
gives 



1 



ai{k) 



ai 2 



(A3) 



with a; the zero-energy scattering length and r; the ef- 
fective range. 

For a square- well potential with depth Vb(< 0) and 
range ro, in the following we scale all lengths {ai,ri) in 
the unit of rg and all momenta (fc,p) of l/rp; for instance, 
a — a/rQ, k = ka^. We define two functions as ji{x) = 
x^^ii{x), ni{x) = x''^^ni{x), modified respectively from 
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the spherical Bessel and Neumann functions. We then 
have 



tan (5; ^ ji-iip)jiik) - ji-iik)jiip) 



(A4) 



with k = y/rnE, p = ^Jm{E — Vq). Note that for / = 0, 
we have j-i ~ —n^^ n^i — jo- The zero-energy scatter- 
ing length and effective range are given by (g = ^J—mVo) 



1 



ai = 



(2Z-l)!!(2; + l)!!jz_i(g)^ 
1 



fi = (2; -1)!!(2; + !)!!{- 



21 



(A5) 

iji-i(g) 



21 - 1 



2l + 3'ji+i{q) 



l(9)^2| 



(A6) 



Appendix B: Scattering under a square- well 
potential with isotropic SO coupling 

Using the continuity properties of ipOj "01 ^-^d their 
first-order derivatives at the potential boundary r = tq, 
we obtain four coupled equations which can be expressed 
in a matrix form. For convenience, we scale all momenta 
in unit of l/ro as A: = kao- 

For scattering state {E = k"^ /m > 0), the matrix equa- 
tion is A{1, t, a, P)^ = with matrix 



q2ji{q2] 

\ 92 jo (52; 



jo (gi ) _ jo (%) - tan 8na (fc^) ^ jo (fc^) 4 

9iji(5i) fc2(ji^fc2) - tanfei^fca)) fci(jipi) 

ji (?L ) _ ji (%) - tan (5ni (fc2) ^ ji (fcj^) 4 

5ijo(gi) fc2(jo(fc2) - tan57io(fc2)) fci(jo(fci) 



tan(5?io(fci) \ 
f tan 5ni(ki)) 
tan i5ni(A;i) 
f tan(5no(fci)) j 



(Bl) 



here we have used jo(a;) = -ji(a;), j((a;) = -f ji(a;) + 

jo (a;) to simplify the equations, k = V m_E, q 

A:2 = A -I- fc, fci = A — fc; 92 = A + V^P^T^, = A — 

The zero determinant |^| =0 gives rise to two solu- 
tions of phase shift 5. When A = 0, these two solutions 
are respectively resulted from two decoupled equations, 
and reproduce the well-known s-wave and p-wave phase 



shifts in free space as given by Ea lA4l 

For bound state, the equations can be obtained 
straightforwardly by transformations (see Ea l43p from 
the equations of scattering state. The binding energy 
E = —K^/m < can be determined from the resulted 
matrix equation \Ai,\ ~ 0. When A = 0, the binding 
energies respectively reproduce the s-wave and p-wave 
results without SO coupling. 
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